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ABSTRACT 


In  this  note  we  consider  the  linearized  equation  for  the  buckling  of  a 
viscoelastic  rod  from  an  undeformed  virgin  state.  We  show  that  this  equation 
does  not  exhibit  buckled  solutions  for  axial  end  thrusts  which  -  after 
application  -  are  held  constant.  Though  this  result  is  apparently  known, 
there  appears  to  be  no  proof  available  in  the  literature. 

We  show  further  that  if  the  load  is  allowed  to  vary  with  time,  then,  in 


contrast  to  elastica  theory,  there  is  an  uncountably  infinite  number  of 
buckled  solutions. 
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SIGNIFICANCE  AND  EXPLANATION 


The  recent  use  of  polymeric  materials  for  structural  purposes  renders 
important  the  careful  study  of  viscoelastic  buckling,  a  phenomenon  quite 
different  from  that  associated  with  elastic  materials,  especially  when  the 
time  scale  of  interest  is  large.  In  this  paper  we  show  that  the  linearized 
equation  for  the  buckling  of  a  viscoelastic  rod  from  its  virgin  state  does  not 
exhibit  buckled  solutions  for  axial  end  thrusts  which  -  after  application  - 
are  held  constant.  We  show  further  that  if  the  load  is  allowed  to  vary  with 
time,  then,  in  contrast  to  elastica  theory,  there  is  an  uncountably  infinite 
number  of  buckled  solutions. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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ON  THE  BUCKLING  OF  LINEAR  VISCOELASTIC  RODS 

11  2 
Morton  E.  Gurtin  ,  Victor  J.  Mizel  ,  and  David  W.  Reynolds 

1.  Basic  equations . 

Consider  a  thin  inextensible  rod,  pinned^  at  the  ends,  and 
acted  on  by  an  axial  end  thrust  P(t)  in  such  a  way  that  its 
center- line  bends  in  a  plane.  For  each  material  point  x  and 
time  t,  let  cp(x,t)  denote  the  angle  between  the  horizontal 
axis  and  the  tangent  to  the  rod  at  x  (Figure  1) ,  and  let 
m(x,t)  designate  the  bending  moment  at  x.  Here,  for  con¬ 
venience,  we  label  material  points  by  their  positions  x  e  [0,L] 
in  the  undeformed,  straight  configuration  the  rod  is  assumed  to 
be  in  prior  to  time  zero. 

We  work  within  the  quasi-static  theory;  thus  balance  of 
4 

moments  has  the  form 

m'  +  P  sin  cp  «  0,  (1) 

where  m1  =  dn/dx.  We  assume  that  the  rod  is  linearly  visco¬ 
elastic  in  the  sense  of  the  constitutive  equation^ 

t 

m(x,t)  =  0cp'(x,t)  +  J  a(t-s)q)  *  (x,s)ds  (2) 

0 

giving  the  bending  moment  as  a  function  of  the  past  history  of 
the  curvature  cp».  Here  0  >  0  is  the  instantaneous  flexural 
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3  In  the  interest  of  brevity,  we  consider  only  pinned  ends.  Our 
results  go  through,  without  change,  for  any  of  the  standard 
boundary  conditions.  (In  this  connection,  cf.  [1].) 
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Cf .,  e.g.,  [2]  ,  §262. 

■’Here  we  use  the  fact  that  cp(x,t)  =  0  for  t  <  0. 
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rigidity,  while  a,  assumed  continuous  and  <0,  measures  internal 
dissipation;  the  unique  solution  G(s),  s  0,  of  the  initial- 
value  problem 

G(s)  =  a(s),  G  (0)  =  0 

is  the  (moment-curvature)  relaxation  function. 

Equations  (1)  and  (2)  yield  the  integro-differential  equa¬ 
tion 

t 

0cp*'(x,t)  +  J  a(t-s)cp"  (x,s)ds  +  P(t)sin  cp(x,t)  =  0, 

0 

or  equivalently,  using  the  standard  notation  for  convolutions, 

0cp"  +  a  *  cp"  +  P  sin  cp  *  0. 

Since  the  ends  of  the  rod  are  pinned,  the  relevant  boundary 
conditions  are 


m(0,t)  =  m(L,t)  =  0 

for  all  t,  or,  by  (2)  and  the  standard  uniqueness  theorem  for 
linear  Volterra  integral  equations, 

cp'(0,t)  =  cp’(L,t)  =  o 

for  all  t.  Further,  as  the  end  of  the  rod  corresponding  to 
the  material  point  x  =  L  can  move  only  horizontally,  we  have 
the  additional  constraint 


-3- 


r 


L 

J  sin  cp(x,t)dx  =  0. 

0 

In  this  note  we  confine  our  attention  to  the  linear  theory 
and  hence  replace  sin  cp  by  cp  in  the  above  equations.  The 
resulting  boundary-value  problem  then  takes  the  form'*' 

0cp"  +  a  *  cp"  +  pep  =  o, 

cp»(0,t)  »  cp'  (L,t)  =  0,  (3) 

L 

J  cp  (x,t)  dx  =  0  . 

0 


Cf.  [3] , [4]  for  a  discussion  of  equations  similar  to  (3)^, 
but  with  P  =  constant.  To  our  knowledge  no  previous  paper 
has  allowed  P(t)  to  depend  on  t. 
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2.  Impossibility  of  buckled  solutions  when  P(t)  =  constant . 

In  this  section  we  show  that  the  problem  (3)  has  only  the 
zero  solution  when 

P(t)  =  P  =  constant  for  t  0. 

With  this  in  mind,  let  0  <  P^  <  < . .  .<  Pn  <••.  and 

^ (x) ,^2 (x) , . . . ,$n(x) , . . .  denote  the  eigenvalues  and  eigen¬ 
vectors  of  the  corresponding  elastic  problem: 

W;  +  pa*n  =  °' 

K{0)  =  K(L)  =  °»  (4) 

L 

J  i()n(x)dx  =  0; 

0 

i.e., 

2  2 

</)n(x)  =  cos  ,  Pn=3£-:|-,  n  =  1,2,...  (5) 

L 

Suppose  that  cp(x,t)  is  a  sufficiently  smooth  solution  of  (3) 
and  define 

L 

cpn(t)  =  J  cp(x,  t)0n(x)dx.  (6) 

0 

Then,  since  cp  and  0n  have  vanishing  spatial  derivatives  at 
the  ends  of  the  rod,  two  integrations  by  parts  in  conjunction 
with  (4)  yield 
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J  ®",6na*  =  J  -  -  tP  V 

0  0 


Thus  multiplying  (3)^  by  and  integrating  with  respect  to 

x  over  the  interval  (0,L)  results  in  the  integral  equation 


<P-Pn)tl’n  -  f  “  *  =  0 


(7) 


for  n  =  1,2.... 

If  P  /  Pn,  (7)  is  a  Volterra  equation  of  the  second  kind, 
which  has  the  unique  solution 


cpn(t)  =  0 


(8) 


for  all  t.  On  the  other  hand,  if  P  =  PR, 

a  *  cpn  ®  o, 

and  hence,  by  Titchmarsh's  theorem1,  (8)  holds  in  this  case  as 
well.  By  (6)  and  (8) ,  cp  is  orthogonal  to  each  function 
of  the  form  (5)  .  This  is  clearly  possible  only  if  cp(x,t) 
is  independent  of  x,  and  the  desired  conclusion,  cp  s  0, 
follows  from  (3)3« 


^Cf.,  e.g.,  [5],  Theorem  152. 
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3.  Existence  of  buckled  solutions  with  P(t)  £  constant . 

There  are  buckled  solutions  for  certain  nonconstant  loads. 
To  see  this  consider  solutions  of  the  form 

cp(x,t)  =  y(t  )^n(x) 


for  t  0.  By  (4)  2  3  this  function  satisfies  the  boundary 
conditions  (3)2  and  the  constraint  condition  (3)3-  Moreover, 
by  (4)-^,  cp  will  satisfy  the  differential  equation  (3)^  pro¬ 
vided  Y  and  P  satisfy  the  integral  equation 

p  i 

[P(t)  -  Pn]Y(t)  -  -f  j  a(t-s)  Y  (s)ds  =  0.  (9) 

0 

This  equation  can  have  many  solutions.  A  simple  example 
is 


Y (t)  =  constant,  P(t)  =  PnQ(t), 


with 


Q(t) 


_  om 

"  6(0)  ' 


For  this  solution  the  angle  cp(x,t)  jumps  to  its  elastic  shape 
at  t  =  0  and  remains  there  for  t  >  0.  The  requisite  axial 
thrust  is  initially  the  elastic  buckling  load  P  ,  but  for 
t  >  0,  P(t)  decreases  in  proportion  to  the  relaxation  function 
G  ( t)  . 
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Solutions  continuous  in  time  at  t  =  0  are  also  possible. 


For  example, 

t 

Y(t)  =  At  (A  =  constant),  P(t)  =  pn  t  J  Q(s)ds> 

0 

forms  a  solution.  (Note  that,  again,  P(0+)  =  pn*) 

More  generally,  (9)  can  be  rewritten  as 

1M  =  V(t)  +  y  *  Q  (t)  . 

Pn  Y(t) 

thus  any  function  y  >  0  on  [0,cd)  yields  a  buckled  solution 

with  P  given  by  (10)  .  it  is  not  difficult  to  show  that 

+  4*  *  + 

P (0  )  exists  when  y (0  )  =  0,  provided  y(0  )  >  0.  Thus  any 

smooth  function  y  on  [0,a>  )  with  y  >  0  on  (0,oo),  y(0)  =  0 

*  +  1  ■ 
and  y (0  )  >0  also  generates  a  buckled  solution. 

Thus,  interestingly,  in  contrast  to  elastica  theory 

there  is  an  uncountably  infinite  number  of  buckled  solutions 

(even  modulo  multiplicative  constants) . 
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^Conditions  on  P(t)  which  insure  the  existence  of  nontrivial 
y (t)  are  given  in  [1] .  This  paper  also  investigates  the 
asymptotic  behavior  of  buckled  solutions . 
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